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ABSTRACT 

We have revised the SWIRE Photometric Redshift Catalogue to take ac- 
count of new optical photometry in several of the SWIRE areas, and in- 
corporating 2MASS and UKIDSS near infrared data. Aperture matching 
is an important issue for combining near infrared and optical data, and we 
have explored a number of methods of doing this. The increased number 
of photometric bands available for the redshift solution results in improve- 
ments both in the rms error and, especially, in the outlier rate. 

We have also found that incorporating the dust torus emission into 
the QSO templates improves the performance for QSO redshift estima- 
tion. Our revised redshift catalogue contains over 1 million extragalactic 
objects, of which 26288 are QSOs. 

Key words: infrared: galaxies - galaxies: evolution - star:formation - 
galaxies: starburst - cosmology: observations 



1 INTRODUCTION 

Rowan-Robinson et al (2008) reported photometric 
redshifts for over 1 million galaxies in the Spitzer 
SWIRE survey (the SWIRE Photometric Redshift 
Catalogue, hereafter SPRC) and gave a detailed re- 
view of earlier work on photometric redshifts. Sub- 
sequently new results have been reported by Wolf et 
al (2008), Brammer et al (2008), Ilbert et al (2009), 
Salvato et al (2009), and a comparison of photomet- 
ric redshift methods has been published by by Hilde- 
brandt et al (2010). 

The SWIRE Survey consisted of 49 sq deg of sky 
surveyed by Spitzer at 3.6, 4.5, 5.8, 8.0, 24.0, 70.0 and 
160.0 fim. The SWIRE areas in which we had opti- 
cal photometry and were able to derive photometric 
redshifts in SPRC were: (1) 8.72 sq deg of ELAIS-Nl 
(ENl), in which we had 5-band (U'g'r'i'Z') photom- 
etry from the Wide Field Survey (WFC), (2) 4.84 
sq deg of ELAIS-N2 (EN2), in which we had 5-band 
(U'g'r'i'Z') photometry from the WFS, (3) 7.53 sq 
deg of the Lockman Hole, in which we had 3-band 



photometry (g'r'i') from the SWIRE photometry pro- 
gramme, with U-band photometry in 1.24 sq deg, (4) 
4.56 sq deg in Chandra Deep Field South (CDFS), 
in which we had 3-band (g'r'i') photometry from 
the SWIRE photometry programme, (5) 6.97 sq deg 
of the XMM Large Scale Structure Survey (XMM- 
LSS), in which we had 5-band (UgriZ) photometry 
from Pierre et al (2007). (6) 1.5 sq deg of ELAIS-Sl, 
in which we have 3-band (B,V,R) photometry from 
Berta et al (2007). In addition within XMM we had 
10-band photometry (ugrizUBVRI) from the VVDS 
programme of McCracken et al (2003), Le Fevre et al 
(2005) (0.79 sq deg), and very deep 5-band photome- 
try (BVRi'z') in 1.12 sq deg of the Subaru XMM Deep 
Survey (SXDS, Furusawa et al 2008). The SWIRE 
data are described in Surace et al (2004) ^. The esti- 
mated redshifts range from to 6 but the bulk are at 
z < 1.5. 



^ and at 
http://irsa.ipac.caltech.edu/data/SPITZER/SWIRE/ 
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The advent of revised INT WFC UgriZ opti- 
cal fluxes for Lockman, ENl and EN2 (Gonzalez- 
Solares et al 2011), the release of CFHT Legacy Sur- 
vey T0005 MegaCam ugriz optical fluxes for XMM 
(littp://tcrapix.iap.fr/cplt/oldSitc/Dcscart/CFHTLS- 
T0005-Release.pdf) and the UKIDSS DR8 release of 
WFCAM JK fluxes for Lockman, ENl and XMM-LSS 
(Lawrence et al. 2007 ^), make it worthwhile revisiting 
the SWIRE photometric redshifts in these areas. 

The Spitzer-solectcd 'data fusiori'catalogucs com- 
piled by Vaccari et al. (2012 in prep) comprises most 
publicly available photometric and spectroscopic data 
such as the above in most fields surveyed by Hcrschcl 
as part of the HcrMES survey (Oliver ct al. 2012). 
TOPCAT was used to merge the SPRC catalogues 
for Lockman, ENl, EN2 and XMM-LSS with the data 
fusion catalogues, with a 1.5" search radius, to get 
the desired photometric data. Not all sources from 
the original SPRC found matches, mostly because the 
data fusion selection requires a source to be detected 
at either 3.6 or 4.5 /um and some 24 //m sources from 
the original SPRC are thus missing from the data fu- 
sion catalogues. Table 1 summarizes the optical and 
near infrared data available in each SWIRE area. 

Rowan-Robinson et al (2008) used optical mag- 
nitudes and Spitzer IRAC 3.6 and 4.5 /xm fluxes to 
estimate photometric redshifts. They reported diffl- 
culty in incorporating 2MASS and UKIDSS J,II,K 
magnitudes into the solution, which they attributed 
to issues of aperture matching. These problems have 
been solved here. The use of these additional bands, 
together with the improved optical photometry, has 
resulted in a reduction in the number of catastrophic 
outliers and improved rms values, when photomet- 
ric and spectroscopic redshifts are compared. Since 
redshifts are now determined from up to 15 bands, 
compared with generally a maximum of 7 previously, 
the new redshifts are more reliable. The new cata- 
logues comprise 1,009,607 redshifts, out of a total of 
1,066,879 in the original SPRC. 



2 APERTURE CORRECTIONS 

Aperture matching between wavelength bands is cru- 
cial to the success of photometric redshifts. For dis- 
tant galaxies photometry in a 2 or 3 arcsec aperture 
will give the integrated light from the whole galaxy. 
For nearby galaxies photometry in the same aperture 
would be dominated by light from the central regions 
of the galaxy and might comprise only a few % of the 
integrated light. In template-fitting methods such as 
we are using here it is natural to try to seek an es- 
timate of the integrated spectral energy distribution 
(SED), so that near and distant galaxies can be fitted 
with the same template. This also has the benefit that 

^ http://surveys.roe.ac.uk/wsa/dr8plus-release.html 
^ http://www.star.bris.ac.uk/~mbt/topcat/ 



derived properties such as luminosity, star-formation 
rate, stellar mass and dust mass have a physical mean- 
ing for the galaxy. 

There are several options for estimating the inte- 
grated light from an extended galaxy in any particular 
waveband. Optical and near infrared catalogues gen- 
erally provide Kron and Petrosian integrated magni- 
tudes. Sextractor (Bertin and Arnouts 1996) provides 
a mag-auto integrated magnitude (similar to Kron 
magnitude). These integrated magnitudes are derived 
essentially by a curve of growth fitted to photometry 
derived in a series of apertures of different sizes. How- 
ever in practice using integrated magnitude estimates 
for each photometric band to derive the integrated 
SED gives poor results for photometric redshifts. This 
is presumably because of the uncertainty introduced 
by the process of estimating the integrated magni- 
tude, primarily because of different contributions of 
sky photon noise. A much more successful option is to 
start from photometry derived in a single small aper- 
ture in each band, and then apply an aperture correc- 
tion derived in a single chosen band to all the bands. 
This is the approach followed by Rowan-Robinson et 
al (2008) and Gonzalez-Solares et al (2011). 

2.1 Optical data 

The optical photometry in SPRC had been derived in 
most areas using Sextractor and we used magnitudes 
measured in a 2" diameter aperture, applying an r- 
band aperture correction 

delmagspijc = r(mag-auto) - r, (1) 

where r is measured in a 2" aperture (WFC 

'aper2'), 

to all bands. All optical magnitudes used in this 

work are AB magnitudes and all J,II,K magnitudes 
are Vega magnitudes (with the exception of the VVDS 
area, where they are AB magnitudes). The quoted 
optical and near infrared magnitudes have been PSF 
(point spread function) aperture corrected, to correct 
for the light lying outside the aperture for a point 
source. Hereafter in this paper we use the term 'aper- 
ture correction' to mean the extended source aperture 
correction. Fig IL shows the optical aperture correc- 
tion, delmagspHc, used in SPRC, versus redshift for 
the Lockman area, essentially a plot of (inverse) angu- 
lar size versus redshift. The aperture correction was 
only applied to the optical magnitudes if it lies in 
the range -0.10 to -5.0. Otherwise it was set to zero. 
Sources with aperture correction > -0.10 mag are con- 
sidered to be point sources within the uncertainties of 
the photometric accuracy and the estimation of the 
integrated magnitude. Sources with aperture correc- 
tion < -5.0 mag. (< 0.01% of sources) are considered 
to be erroneous determinations of integrated magni- 
tude. 
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Table 1. Explanation of magnitude fields 



field 


aml-5 


am21,22,23,25,26 


am6-8 


ENl 


rev. WFC UgriZ 


SDSS ugriz 


2MASS/UKIDSS JHK 


EN2 


rev. WFC UgriZ 


SDSS ugriz 


2MASS JHK 


Lock 


rev. SWIRE Ugri 


SDSS ugriz 


2MASS/UKIDSS JHK 


VVDS 


VVDS ugri 


UBVRI 


VVDS JK 


SXDS 


Subajru 


Mcgacam ugriz 




XMM 


SWIRE UgriZ 


Mcgacam ugriz 


UKIDSS JHK 


CDFS 


SWIRE gri 






SI 


SWIRE B,V,R 







For SDSS optical data, available for the Lock- 
man, ENl and EN2 areas, the 'model' magni- 
tude provides a well-calibrated integrated magnitude. 
Fig IR shows delmagspi?c versus r(SDSS,model)- 
r(WFC,mag-auto). Quasars almost all have no aper- 
ture correction, as expected. For galaxies there 
is no sign of any correlation of r(SDSS,model)- 
r(WFC,mag-auto) with delmagspi?c, so the SDSS 
model magnitude is estimating approximately the 
same total magnitude as the WFC mag-auto. 
However there is quite a wide dispersion in 
r(SDSS,model)-r(WFC,mag-auto), sufficient to harm 
photometric redshift estimation (see SED plots be- 
low, Fig 14). This dispersion arises from the uncer- 
tainties in the two estimates of the integrated mag- 
nitude. We have therefore added r(WFC,mag-auto)- 
r(SDSS,model) to the SDSS model magnitudes to nor- 
malize the two sets of aperture corrections. The rea- 
son for making the WFC r-magnitude the preferred 
choice is that there are far more SWIRE sources with 
WFC data than with SDSS data. Also the WFC 
goes one magnitude deeper than SDSS and there- 
fore provides more accurate magnitudes for fainter 
SDSS sources. For SDSS data we investigated using 
the point-source (PS) magnitudes, with aperture cor- 
rection r(SDSS,Petr)-r(SDSS,PS), but this gave infe- 
rior results to using the SDSS model magnitudes. This 
is understandable because the model magnitudes use 
an aperture derived in the r-band and then measure 
the flux in that aperture in all the bands, ensuring 
consistent galaxy colours. 

We investigated various other options for aper- 
ture corrections in the optical. The revised WFC pho- 
tometry of Gonzalez- Solares et al (2011), available in 
the fusion catalogues, includes the Petrosian magni- 
tude in each band, so we define 

delmagwFC = r(WFC,Petr)-r. (2) 

This is well correlated with delmagspi?c (Fig 

2L), but with significant scatter. 

Using delmagvi/FC instead of dehnngs prc led to 
slightly worse phot-z results, so we decided to stick 
with delmagsPHC, which was derived from the Sex- 
tractor r-band mag-auto, for WFC data (while using 
the revised WFC magnitudes supplied in the fusion 



catalogue) . However where a mag-auto estimate is not 
available, we have used delmagwFC. 

2.2 Near infrcired data 

For 2MASS data an option for the aperture correction 
is to use the K-iso magnitude if available, and the K- 
PS magnitude otherwise, with aperture correction 

delmag2MASs=K(2MASS,iso)-K(2MASS,PS). 

(3) 

For UKIDSS data the natural aperture correc- 
tion to consider is 

delmagf/ if/DSS 
=K(UKIDSS,Petr)-K(UKlDSS,aper3), (4) 

applied to the aperS (1.0") J,K magnitudes. 

Both delmag2MASS and debnaguKiDSS are quite 
well correlated with delmagspHc (Fig 2R shows the 
correlation of delmag(//f/DS5 with delmagspRc) and 
this suggests the idea of using k*dclmagspflc as the 
near infrared aperture correction, with k to be deter- 
mined for each survey. The direct use of delmag2MASs 
and delmagr/KiDSS for 2MASS and UKIDSS magni- 
tudes, respectively, resulted in a worse phot-z solu- 
tion, so the use of k*delmagspiic was explored in 
some detail. 

Fig 3 shows delmagsp jic versus (z-J), with no 
aperture correction applied to J. There is a very 
strong correlation. Fig 4L shows delmagsppc ver- 
sus (z-J) for 2MASS data, with an aperture correc- 
tion 0.8*delmagspHc applied to J. Fig 4R shows 
delmagspijc versus (z-J) for UKIDSS data, with an 
aperture correction l.l*delmagsppc applied to J. 
Only sources with spectroscopic redshift <0.3 have 
been included to minimize the possibility of evolution 
of colour with angular size. We can see that these 
corrections work well in removing the correlation of 
colour with delmagsPHC. 

2.3 IRAC data 

In SPRC we used Kron fluxes if the 3.6 /um size was 
greater than a specified threshold (area > 200 pix- 
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Figure 1. LH: SPRC aperture correction versus Zsp for Lockman. RH: SPRC aperture correction versus r(SDSS,model)- 
r(WFC, mag-auto). Red symbols: galaxies, blue symbols: QSOs. 





delmagg, 



-2 
delmagg, 



Figure 2. LH: Plot of Ae\vcia.gw FC (eqn 2) versus delmagsp^c the SPRC aperture correction, in Lockman. The straight 
line has slope 1. RH: Plot of delmagS (eqn 3) versus delmaggpjjc the SPRC aperture correction. The straight line has 
slope 0.8. 



els) , 'IRAC-aper2' fluxes (measured in a 3.8" aper- 
ture) otherwise. For sources smaller than this thresh- 
old, the Kron and aper2 magnitudes agreed well. 
This is also the approach adopted here. Fig 5 shows 
delmagspjjc versus (K-am(3.6/im)) for 2MASS (L) 
and UKIDSS (R) data, where am(3.6/im) is the AB 
magnitude derived from the Spitzer 3.6/im flux. There 
are some residual issues for the UKIDSS-IRAC com- 
parison. This does harm the photometric redshift es- 
timates for some objects, resulting in a larger re- 
duced for the fit. SED fits for selected galaxies 
with > 5 with known spectroscopic redshift shows 
that the IRAC fluxes gives good consistency with the 



aperture corrected optical magnitudes. The problem 
lies with overestimation of the J,H, and especially K 
brightness for extended galaxies. Changing the IRAC 
aperture correction to k*delmagspijc improved the 
appearance of Fig 6R enormously, but led to worse 
phot-z results. The main problem was an offset of the 
log^iil+Zph) points relative to /o(7io(l+2sp) by ~0.01. 
This was not fixed by the process of in-band correc- 
tion factors (Ilbert et al 2006), and it would probably 
be necessary to revise the templates in the near in- 
frared to achieve convergence. 
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2MASS 




UKIBSS 




z-J 



Figure 4. LH: SPRC aperture correction versus corrected (k=0.8) z-J for 2MASS. RH: SPRC aperture correction versus 
corrected (k=l.l) z-J for UKIDSS. Only sources with spectroscopic redshift <0.3 have been included to minimize the 
possibility of evolution of colour with angular size. 




-2 -1 1 -2-10 1 

K-am(3.6) K-am(3.6) 



Figure 5. LH: SPRC aperture correction versus corrected (k=0.7) K-am(3.6) for 2MASS. RH: SPRC aperture correction 
versus corrected (k=l.l) K-am(3.6) for UKIDSS. Only sources with spectroscopic redshift <0.3 have been included to 
minimize the possibility of evolution of colour with angular size. 



3 PHOTOMETRIC REDSHIFTS 

The photometric redshift method is as described in 
SPRC. It is a two-pass template method based on 6 
galaxy ( 11 in the second pass) and 3 AGN templates. 
The galaxy templates were first generated empirically, 
then modelled with a full star-formation history so 
that star-formation rates and stellar masses can be 
derived. After the first pass the infrared excess rela- 
tive to the starlight fit is modelled with four infrared 
SED types (cirrus, M82 starburst, Arp220 starburst 
and AGN dust torus). Dust extinction with Av rang- 



ing up to 1.0 is permitted. For QSOs SMC dust prop- 
erties are assumed. 

For galaxies and QSOs, JHK magnitudes, and 3.6 
and 4.5 ^im fiuxes were used in the first pass through 
the data provided there was no evidence of a strong 
AGN dust torus component, determined by condition 
S(5.8) > 1.2 S(3.6). For galaxies these near infrared 
data were not used in the second pass fit if the dust 
torus component was dominant at 8 /im from the ir 
template fitting. For the second pass for QSOs, a new 
approach was developed, described in section 4 below. 
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z-J 



Figure 3. SPRC aperture correction versus z-J, with no 
aperture correction applied to J. 

to allow use of the near infrared data even though a 
dust torus component is present. 

In the photometric redshift solution we used the 
new SDSS 'model' magnitudes, and the revised WFC 
magnitudes, but we retained the WFC star/galaxy 
classification in each band (because this information 
is available for a much higher proportion of sources) 
and the optical aperture correction used in SPRC. 
We used 2MASS JHK (PSC), where available, and 
UKIDSS JK (aperS) if not. 

We used galaxies with known spectroscopic red- 
shifts to determine in-band correction factors, follow- 
ing Ilbert et al (2006). These are shown in Table 2. 
Spectroscopic redshifts came from a wide range of ref- 
erences (details are given in the on-line catalogue) and 
are mostly in the range 0-1.2 for galaxies. 

Fig 6 shows a comparison of the SDSS logioil -f 
Zphot) with logio{l + Zspect) for the SWIRE Lockman 
sample (LH). and the same plot for the present sample 
(RH). The SDSS performance is better at z < 0.3, but 
our approach works better at z > 0.5. To improve our 
photometric redshifts at z < 0.3 it would be necessary 
to refine our optical templates using the new photo- 
metric data, and to increase the number of templates. 
The latter would be likely to worsen the performance 
at higher redshift through increased aliasing. 

Fig 7L shows the same comparison for the 
SWIRE Photometric Redshift Catalogue, restricted 
to r < 23.5, reduced < 3 and at least 5 photo- 
metric bands in the solution. Fig 7R shows a similar 
plot for the revised SPRC, with the requirement that 
K be selected, and for a minimum of 7 photomet- 
ric bands. We can see that inclusion of near infrared 
(JHK) data in the solution has reduced the number 
of catastrophic outliers. 

Fig 8L shows the same comparison for the re- 
vised catalogue in XMM, with at least 9 bands in 



the solution in the solution. Fig 8R shows the same 
comparison for photometric redshifts from the LeP- 
hare method (Ilbert et al 2006). The latter results are 
better for z < 0.5, but worse for z > 0.5. Again this 
is probably the consequence of the larger number of 
templates used by Ilbert et al. 



4 QSOS 

The SPRC approach required that an object be 
flagged as stellar to consider a QSO optical template. 
With the SPRC stellar flag some QSOs get missed 
(and end up with the wrong redshift) as a result of 
this condition. As discussed in SPRC it is not possi- 
ble to allow a QSO template option for all galaxies, 
since far too many galaxies end up with mistakenly 
high redshifts. However, while still requiring a stellar 
flag for an object to select a QSO template, we have 
allowed the SDSS stellar flag to override the WFC 
flag where they disagree and this allows a few more 
quasars through. 

Salvato et al (2009) have demonstrated excellent 
performance for 1032 QSOs and AGN in the COS- 
MOS field, using 30 photometric bands, including 12 
narrow-band filters. They introduce two innovations: 
firstly they track the variabilty of quasars and ap- 
ply an appropriate correction to the photometry. Sec- 
ondly they use templates that include a range of con- 
tributions from AGN dust tori. 

While we do not have the information to track 
QSO variability in our photometry, we have explored 
the idea of adding a range of AGN dust tori strengths 
to our QSO templates during the second pass, and 
then using the 1.25-8 jj, m data in the redshift solu- 
tion. The amplitude of dust tori added corresponded 
to Ltor/Lopt = 0, 0.2, 0.4, 0.6, 0.8, 1.0. The perfor- 
mance for QSOs is much improved (Fig 9), especially 
in the reduction of catastrophic outliers. For at least 
11 photometric bands, reduced < 3, r < 21.5, we 
find an rms of 9.3% and an outlier rate 9.3%. For com- 
parison, Salvato et al (2009), with 30 bands, achieved 
an rms of 1.2% for QSOs with / < 22.5, and an 
outlier rate of 6.3%. Their greatly improved rms can 
be attributed to the correction for variability and to 
the use of 30 photometric bands, including 12 narrow 
band filters. But our outlier performance is almost as 
good, despite less than half the number of photomet- 
ric bands. 



5 REASONS FOR POOR 

Although the photometric redshift estimates are good 
for sources where the reduced < 3, the estimates 
get slightly worse as increases and some sources 
have very large x^- We have investigated the reasons 
for poor x^ through colour-colour diagrams and SED 
plots. The main reasons found for poor are found 
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Table 2. Correction factors for fluxes in each band, by area 



area 


U' 


g' 


r' 


i' 


Z' 


J 


H 


K 


3.6 /jm 


4.5 fim 


EN1,EN2 


f.354 


f.l35 


f.l04 


1.059 


0.650 


0.976 


0.900 


0.888 


1.071 


1.104 


Lockman 


f.360 


f.062 


0.982 


0.922 




0.923 


0.904 


0.899 


1.036 


1.066 


VVDS 


f.046 


0.969 


0.928 


0.960 


0.972 


1.173 




1.029 


1.085 


1.015 


XMM-LSS 


f.308 


f.206 


f.056 


0.925 


1.011 








1.085 


1.086 


CDFS 


f.046 


0.969 


0.916 


0.991 


0.972 








1.238 


1.037 




Figure 6. LH: SDSS photometric redshifts versus spectroscopic redshifts for Lockman SWIRE sample. RH: photometric 
redshifts from present work versus spectroscopic redshifts for Lockman SWIRE sample. 



Ozl 234 Ozl 234 




Figure 7. LH: Photometric redshifts from SWIRE Photometric Redshift Catalogue (Rowan-Robinson et al 2008). for 
galaxies with r < 23.5, at least 6 photometric bands, reduced < 3. Red symbols: galaxies. RH: Same for revised SPRC, 
with at least 11 photometric bands. 
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Figure 8. L: Photometric redsliifts in XMM-LSS using new fusion catalogue: ugriz from SDSS, revised ugri from WFC, 
JHK from 2MASS, JK from UKIDSS. R: Photometric redshifts from LePhare method (Ilbert et al 2009). 




Figure 9. LH: Photometric redshifts from SWIRE Photometric Redshift Catalogue (Rowan-Robinson et al 2008). for 
quasars with r < 23.5, at least 7 photometric bands, reduced < 3. Blue symbols: QSOs. RH: Same for revised SPRC, 
with at least 12 photometric bands. 



to be (i) contamination by stars, (ii) poor photom- 
etry, (ill) QSOs misclassified as galaxies, (iv) QSOs 
for which the extinction is higher than the maximum 
assumed Av = 1. 

Stars can be recognised in a plot of the ratio of 
3.6 /^m to r-band flux versus r-i ( Fig lOL, see Rowan- 
Robinson et al 2005). A clear stellar sequence can be 
seen. We have removed stars from the catalogue by 
excluding 17326 sources with stellar flags which lie 
below the solid line in Fug lOL, and which had re- 
duced > 7. Figure lOR shows the ^nd pho- 
tometric redshift generated for these candidate stars 
by the photometric code. Stars can generate spectac- 



ularly bad x^ a-nd tend to alias with specific galaxy 
and QSO redshifts. We believe that few of the sources 
with <7 are stars. 

Figure IIL shows (u-g) versus (g-r) (WFC pho- 
tometry) for QSOs, with blue symbols for sources 
with x^ < 3 and red symbols for x^ > 7, after removal 
of stars using the 3.6/^m-(r-i) constraint. The bro- 
ken locus is for blackbodies with temperature ranging 
from 3000 to 15000, the solid locus shows the effect 
of extinction using SMC dust properties. The > 1 
objects show a much wider range of colours than those 
with x^ < 3. Examination of SED plots shows that 
most of the x^ > 7 objects are stars with slightly 
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j. :K; 

■ i.'l- 



0,4 0,6 0.8 

Figure 10. LH: Ratio of 3.6 /xm to r-band flux versus r-i for whole catalogue, for sources with > 20 and at least 12 
photometric bands. A clear stellar sequence is seen. Stars have been excluded from the catalogue by removing objects 
flagged as point sources lying below the solid line. RH: Reduced versus estimated photometric redshift for sources 
classifled as stars. Red symbols: sources which the photometric redshift code identified as galaxies; blue symbols: QSOs. 



-2 




5 (Av,sMc) 



errl<3 




-2 



Figure 11. LH: (u-g) versus (g-r) (WFC photometry) for QSOs in ENl field with at least 11 photometric bands, with blue 
symbols for sources with < 3 and red symbols for > 7, after removal of stars using the 3.6/tm-(r-i) constraint. The 
broken locus is for blackbodies with temperature ranging from 3000 to 15000, the solid locus shows the effect of extinction 
(Ay=0-5) assuming SMC dust properties. RH: Reduced versus estimated photometric redshift for sources classified as 
stars. Red symbols: sources which the photometric redshift code identified as galaxies; blue symbols: QSOs. 



higher 3.6 fim fluxes than the stellar sequence in Fig 
lOL. A few are QSOs with higher extinction than the 
maximum Ay = 1 we have used. Others have high 
simply because of poor photometry or SED shapes 
differing from our templates. 



Figure IIR shows (u-g) 
ies with SDSS photometry, 



versus (g-r), for galax- 
with difi^erent symbols 



for X < 3 and > 7. About 10% of 'galaxies' with 
X'^ > 7 achieved much improved fits with QSO tem- 



plates (some of these did not have the stellar flag that 
was a prerequisite for consideration as a QSO). Some 
are stars with slightly higher 3.6 ^im fluxes than the 
stellar sequence in Fig lOL. For the rest the problem 
appears to be poor photometry (and underestimated 
photometric errors). 

Figure 12L shows (i-z) versus (r-i) (WFC pho- 
tometry) for galaxies in ENl with at least 11 photo- 
metric bands, with different symbols for x^ < 3 and 
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Figure 12. L; (i-z) versus (r-i), WFC photometry, for galaxies witli at least 11 photometric bands and > 7 (red), < 3 
(blue). R: Histogram of reduced for ENl area, after exclusion of stars. 



> 7. The galaxies with higher show a wider spread 
in colours, suggesting the main issue is photometric 
problems. Figure 12R shows a histogram of the re- 
duced for sources in the ENl area, after removal 
of stars. 



6 REVISED SWIRE PHOTOMETRIC 
REDSHIFT CATALOGUE 

Our revised SPRC contains 818555 objects, 204421 in 
ENl, 116195 in EN2, 217461 in Lockman, and 280478 
in XMM-LSS, compared with 875353 in the same ar- 
eas in the original SPRC. 3.6% of SPRC sources did 
not find a match in the fusion catalogues, mainly 
because the latter omitted 24 ^m only sources. A 
further 1% failed to achieve a redshift solution, ei- 
ther because there were less than two valid photo- 
metric bands or because the reduced > 100. The 
SWIRE redshifts in CDFS and SI have not been re- 
vised because we have no new photometric informa- 
tion in these areas. These two areas bring the to- 
tal number of redshifts in the revised catalogue to 
1009607, with readmeSWIRErev in same directory). 
Our new catalogue delivers photometric redshifts for 
26288 quasars. 

Fig 13L shows a histogram of the number of pho- 
tometric bands for the new catalogue, compared with 
that for the SPRC. Fig 13R shows the redshift dis- 
tribution of galaxies and quasars in the new cata- 
logue. Although most of the galaxy redshifts are at 
z < 1.5, there is a broad tail to higher redshift. Our 
catalogue is a very significant resource for galaxies at 

4 

http: / /astro. ic.ac.uk/ '^mrr/swircphotzcat/zcatrevl2ff3.dat.gz 



high redshift, with 106537 sources with z > 2 and 
15226 sources with z > 4. 



7 SEDS OF OUTLIERS 

To investigate outliers, we have plotted in Fig 14 the 
SEDs of 8 outliers from the Zphot v. Zapec comparisons, 
those with < 3, and with more than 12 photomet- 
ric bands in ENl, EN2 and Lockman, and with more 
than 9 bands in XMM-LSS. SDSS photometry (af- 
ter applying the aperture correction of section 2.1) is 
shown in red. Details of the sources are given in Ta- 
ble 3. Two objects have aliases at the spectroscopic 
redshift, one appears to need an Sab template with 
Av =0.35 {Av is set to zero for the Sab template in 
our code), and 4 could be errors in the spectroscopic 
redshift (two were based on a single line). 



8 DISCUSSION 

With the changes discussed in the previous sections, 
we now see improvements in the fraction of catas- 
trophic outliers in the photometric redshift estimates. 
The rms errors for galaxies are very similar to those 
achieved in SPRC, whole those for quasars are signif- 
icantly improved. 

For galaxies with reduced x^ < 3, r < 23.5, the 
rms in {zphot — Zspec)/{1 + Zspec), after rejection of 
outliers with values discrepant by 15% or more, is 3,7, 
3.4 %, for no. of bands 10, 14 respectively. The cor- 
responding percentages of outliers are 1.2 and 0.2 %, 
respectively. Fig 15L shows how the percentage of out- 
liers for galaxies with |/opio(l -I- Zphot)/ii + Zspect)\ > 
0.06, vary with the number of photometric bands. Fig 
15R shows the percentage rms {a^/il + z)) versus 
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Figure 13. L: Histogram of number of bands used in estimating photometric redshifts in revised catalogue (solid line) 
compared with SPRC (broken line). R: Histogram of redshifts in revised SWIRE Photometric Redshift Catalogue (red: 
galaxies, blue: QSOs, x 10). 



Table 3. Zpf^ V. Zsp outliers whose SEDs are plotted in Fig 10, from bottom 

no. of bands notes 



RA 


dec 


Zspec 


ref. 


^phot 


no. 


161.30797 


58.74822 


0.2470 


1 


0.52 


12 


161.72533 


58.99032 


0.1890 


1 


0.55 


12 


161.74496 


58.94529 


0.8220 


1 


0.25 


12 


161.94997 


59.29018 


0.5600 


2 


0.17 


13 


36.40352 


-4.39069 


0.1307 


3 


0.85 


9 


36.47719 


-4.36419 


0.8277 


3 


1.24 


9 


248.47900 


41.48283 


0.1592 


4 


0.47 


14 


1 Owen et al 2009 












2 Rowan-Robinson 


et al 2008 










3 Lc Fevre et al 2005 












4 Rowan-Robinson 


et al 2004 











alias 
alias 

zspec error ? 

zspec error ? 

zspec error (single line) ? 

zspec error (single line) ? 

need Sab with Av=0.35 



number of photometric bands. The additional photo- 
metric bands provided by the JHK data from 2MASS 
and UKIDSS now have a clear beneficial effect, espe- 
cially on the fraction of catastrophic outliers. Our re- 
sults can be compared directly with the corresponding 
results for SPRC, and with new results from EAZY 
(Brammer et al 2008), GOODS-S (Dahlen et al 2010), 
and COSMOS30 (Ilbert et al 2009). While these com- 
parisons can be affected by differences in the survey 
depths, the main differences between these surveys is 
the number of photometric bands. Our outlier perfor- 
mance is consistently better than other methods. 



Our new SWIRE redshift catalogue should be 
useful for improved studies of the infrared extragalac- 
tic population and for the Herschel surveys carried out 
in all of the SWIRE fields by the Hermes consortium. 
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Figure 14. SEDs of outliers in Zpf^ v. Zsp plot, using the spectroscopic redshift, with details of sources given in Table 1. 
Outliers with more than 12 photometric bands in ENl, EN2, Lockman, and more than 9 bands in XMM-LSS, are included. 
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Figure 15. L: Percentage of outliers versus number of photometric bands for SPRC (black loci) and for fusion catalogue 
(red loci), for sources with reduced < 3, r < 23.5. R: Percentage rms (^^/(l + z)) for SPRC (black loci) and for fusion 
catalogue (red loci), for sources with reduced < 3, r < 23.5. 
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